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Scaling expressions for the free energy are derived, using the Luttinger-Ward (LW) functional 
approach in the Eliashberg framework, for two different models of quantum critical point (QCP). 
First, we consider the spin-density-wave (SDW) model for which the effective theory is the Hertz- 
Moriya-Millis (HMM) theory, describing the interaction between itinerant electrons and collective 
spin fluctuations. The dynamic of the latter arc described by a dynamical exponent z depending 
on the nature of the transition. Second, we consider the Kondo breakdown model for QCP's, one 
possible scenario for heavy-fermion quantum transitions, for which the effective theory is given by a 
gauge theory in terms of conduction electrons, spinous for localized spins, holons for hybridization 
fluctuations, and gauge bosons for collective spin excitations. For both models, we construct the 
thermodynamic potential, in the whole phase diagram, including all kinds of self-energy corrections 
in a self-consistent way, at the one loop level. We show how Eliashberg framework emerges at this 
level and use the resulting Eliashberg equations to simplify the LW expression for free energy . it is 
found that collective boson excitations play a central role. The scaling expression for the singular 
part of the free energy near the Kondo breakdown QCP is characterized by two length scales : one 
is the correlation length for hybridization fluctuations, and the other is that for gauge fluctuations, 
analogous to the penetration depth in superconductors. 

PACS numbers: 71.10.Hf, 71.30.-t-h, 71.10.-w, 71.10.Fd 



I. INTRODUCTION 

Fluctuation corrections are an essential ingredient near 

quantum critical points (QCPs). It may be relatively 
easy to incorporate quantum corrections in the weak cou- 
pling approach, the so called Hertz-Moriya-Millis (HMM) 
theoretical framework"'^. However, it becomes more com- 
plicated to include quantum fluctuations in the strong 
coupling approach such as the gauge theoretical frame- 
work sometimes proposed to describe strongly correlated 
electrons like doped Mott insulators^ and some heavy- 
fermion QCPs'^. For these models, it is believed that 
strong correlations fractionalize electrons into some ex- 
otic elementary excitations carrying fractional quantum 
numbers of electrons, and quantum fluctuations of such 
enhanced degrees of freedom appear to be complicated. 
It is challenging to develop a systematic approach to 
introduce, self-consistently, physically essential fluctua- 
tions into the thermodynamic potential near a QCP. 

Effects of quantum corrections on the thermodynamic 
potential can be incorporated systematically using the 
Luttinger-Ward (LW) functional approach^"^, where the 
grand potential is written in terms of dynamic quantities, 
such as the fully dressed Green's function G[E] and the 
self-energy S, through the relation^ 

f2[S] = T STr[ln{-G-i[E]} + EGp]] -hFjcp]}, (1) 

where STr[A]=Tr[AB]-Tr[A_F] is the supertrace over Mat- 
subara frequencies, internal quantum numbers of the 
bosonic (B) and fermionic (F) components of A. The 



quantity f|g[S]| is the so-called LW functional, deter- 
mined purely by the interaction potential and given by 
the sum of all closed-loop two-particle irreducible skele- 
ton diagrams in the perturbation theory approach. Vari- 
ation of the LW functional Y with respect to G generates 
the self-energy 
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where Gq is the non-interacting Green's function. 

The thermodynamic potential is stationary with re- 
spect to changes of the self-energy, i.e. it satisfies the 
saddle-point condition 



0. 



E=S, 



(3) 



An important issue in the perturbation approach of the 

LW functional is to find an explicit functional dependence 
for y^. This is generally unknown and it is not always 
possible to sum the skeleton expansion into a closed form 
for Y. The problem with strongly correlated systems is 
even worse because the convergence of the skeleton ex- 
pansion is not guaranteed. It was demonstrated that the 
LW functional can be written as a closed form in the 
Eliashberg framework®'^*^, where the Eliashberg approx- 
imation allows to handle quantum corrections in a self- 
consistent way, at the one-loop level. The Eliashberg the- 
ory turns out to be justified for 2; > 1 quantum criticality, 
where z is the dynamical exponent, using a large- A'^ ex- 
pansion supporting the Migdal theorem^^. Here N is the 



2 



number of fermion flavors. In particular, the Eliashberg 
theory was axgued to be a minimal framework well work- 
ing near a QCP^*^. 

In this study, wc derive a LW expression of the free 
energy for two models of itinerant QCPs : the spin- 
density-wave (SDW) model and the Kondo breakdown 
model. This will allow us to describe thermodynamics 
near these QCPs starting from a microscopic model and 
incorporating self-consistently the effect of quantum fluc- 
tuations. 

The plan is as follows. In section II, we present the 

Spin-fermion (SF) model and derive in a systematic way 
the LW functional and show how Eliashberg equations 
for self-energies are derived. The expression of the free 
energy is simplified using Eliashberg equations and the 
scaling expression of its singular part is deduced. Sec- 
tion III is devoted to the Kondo Breakdown (KB) model. 
A particular care is taken to describe the Higgs part of 
the phase diagram. The effect of condensation is incor- 
porated into a zero-order theory before considering a cu- 
mulant expansion in the fluctuations interaction. For this 
gauge theory, there are additional collective excitations, 
which results in the presenc;e of two length scales in the 
scaling expression of the free energy. These two scales are 
related through the Anderson-Higgs mechanism. Section 
IV summarizes and discusses our main results. In partic- 
ular, theoretical structure differences between the HMM 
theory and the gauge theory of the Kondo breakdown 
QCP are emphasized. Technical details are presented in 
the appendices. 



II. REVIEW OF THE LUTTINGER-WARD 
FUNCTIONAL APPROACH IN THE 
ELIASHBERG FRAMEWORK OF THE 
SPIN-FERMION MODEL 

The standard model of quantum criticality in a metal- 
lic system is the HMM theory. In this model, a dynamical 
exponent z, relating the variation of the energy with the 
momentum w ~ g^, characterizes the dynamics of col- 
lective excitations near the QCP. In particular, z = 3 
describes the ferromagnetic QCP while z = 2 describes 
the antiferromagnetic one. It is valuable to review the 
construction of the LW functional in the HMM theoret- 
ical framework, discussed in the past", although several 
heavy fermion compounds have been shown not to follow 
the ^; = 2 HMM theoryi^-so. 
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where we used the "relativistic" notation for energy- 
momentum k = (k,icj), q = (q, ifi) and the sum ex- 
pression is defined as 
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In (4) , tp^rk is the fermionic annihilation field for an 
electron with energy-momentum vector k and spin a, Sg 
is a bosonic field describing spin-fiuctuations near a mo- 
mentum Q and g is the coupling constant measuring the 
strength of the interaction between fermionic and bosonic 
excitations. The last term in (4) stands for higher order 
terms in S. These are shown to be irrelevant for d > 2 
and marginal for d = 2 and can therefore be neglected^^. 

In the absence of the interaction, fermionic and bosonic 
excitations axe described by the bare electron Green's 
function Go(fc) and the bare spin susceptibility Xo(9) re- 
spectively 



Go{k) = T 

xo{q) 
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A. Spin-fermion model 



zq is the quasiparticle renormalization factor given by 
the Fermi liquid theory, and the electron dispersion is 
linearized with a Fermi velocity vp and Fermi momen- 
tum \lf- The electron's chemical potential can incor- 
porate effects of the condensed part of the bosonic field. 
The bare spin susceptibility is the usual Ornstein-Zernike 
form where is the bare correlation length of spins and 
Xo^o is the static susceptibility. Vg is the bare spin veloc- 
ity. 

The spin-fermion model Eq.(4) is an effective low- 
energy model that can be derived from the Hubbard-like 
model in the weak coupling approximation^^'^^'^^; high- 
energy fermions, with energy above A, are integrated out 
to generate collective bosonic modes that mediate the 
interaction between fermions at energies smaller than A. 
Dynamics of the low-energy fermions and the collective 
spin excitations are then described by Eq. (5) . 



We start from the so called spin-fermion model B. Eliashberg theory 

(gp) 10,21,22 £qj. ^Y,e SDW transition 

The Eliashberg framework allows the evaluation of the 
— T'S^ip^ [—GQ^{k)]ipcrk self-energies E and 11, for electrons and spin fiuctuations 

^ respectively, self-consistently assuming we can neglect 
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the momentum dependence of the E and vertex correc- 
tions. An extensive review of this technique is given in^''. 
We recall here the spirit of this technique and main re- 
sults. 

The Ehashberg procedure relies on three steps: 

• neglect both the vertex corrections and the momen- 
tum dependence of the fermionic self-energy : 

E(k,ia;„) = S(ztj„) 
Ag = 

• Use Dysons' equations 

X"^(q,ifi„) = Xo^{(l,i^n) - n(q,ifi„), (6) 

to evaluate self-consistently the self-energies repre- 
sented diagrammatically in Fig. 1(a) and (c), where 
the propagators are fully dressed according to (6). 

• Check a posteriori that the neglected momentum 
dependence of the fermionic self-energy and vertex 
corrections are indeed small. 
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FIG. 1: (a) and (b) are the first order contribution to the 
fermionic self-energy and (c) is the polarization bubble, where 
cr € [i,N] is the spin index. The propagators of the fermion 
(straight line) and the spin fluctuations boson (wavy line) are 
fully dressed, (b) is a static and uniform part in the self- 
energy and can thus be considered as a renormalization of 
the electron chemical potential. 

The bosonic self-energy is found to be 

n(q,za.) = 7^, (7) 

where 7 — g^X^kp ((jivy) and z is the dynamical ex- 
ponent. This result is customary for problems where 
fermions interact with their own collective modes. The 
latter are damped whenever they lie inside the particle- 
hole continuum of the Fermi liquid^. Such a Landau- 
damped term is larger than the regular ©(fi^) term in 



the bare spin susceptibility Eq.(5) and fully determines 
the collective spin dynamics. This causes feedback effects 
on the self-energy correction of electrons, giving rise to 
non- Fermi liquid physics near the QCP"'^°. 

The model (4) can be extended by introducing iV 7^ 1 
identical fermionic species. A channel index v ^ [1,A^] 
is then added to the fermionic operators 1/' in (4) and 
g g/^fN to ensure a well-defined large N limit. 

It has been shown that the Eliashberg approximation 
becomes exact in the limit N 00^'^-'^. Indeed, it is 
shown that both vertex corrections and the momentum- 
dependent corrections to the fermionic self-energy turn 
out to scale as 1 /N and vanish in the limit N ^ 00. This 
limit shares some similarity with the Migdal limit for the 
electron-phonon problem : at large N, the damping in- 
troduced in (7) scales as N and the collective excitations 
become slow. Then, the smallness in 1/iV compares to 
the smallness in m/M where m is the electron's mass and 
M is the ion's mass. 

In the following, we will show how the Eliashberg 
framework emerges from the LW approach to the Spin- 
Fermion model. 



C. Luttinger-Ward functional for the spin-fermion 
model 

As said in the introduction, the LW functional is, di- 
agrammatically, the sum of all closed-loop two-particle 
irreducible skeleton diagrams^. These can be ordered in 
a 1/N expansion as in Fig. 2 
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FIG. 2: Leading skeleton diagrams participating to the LW 
functional Y for the Spin-Fermion model with dependence on 
Fermion (straight line) and boson (wavy line) prop- 
agators are fully dressed, o £ [i,N] is the spin index and 
V £ [1, N] is the channel index. The first diagram contains 
one fermionic loop carrying spin and channel quantum num- 
bers and one pair of vertices, each of order 0{l/y/N), so that 
it is of order 0{N^/N) = 0{N). The second dia gram involves 
one fermionic loop and two pairs of vertices so that it is of 
order 0{N'^/N'^) = 0(1). Bracketed terms are dropped in 
the large N limit. 

Considering the general expression of the LW expres- 
sion [Eq.(l)], and taking into account only the lead- 
ing 0{N) contribution to Y, shown in Fig. 2, one can 



write down the free energy in terms of electron and spin- 
fluctuation Green's functions and self-energies 



Flw[^M = -iVT^[ln{-G-i(fc)} + E(fc)G(A:) 

k 

+ T^[ln{x-i(g)}+n(g)x(<z) 



-I- 3iV5'r2^G(fc)x(g)G(fc + g), (8) 

k,q 

where G{k) and S(fc) are the fully renormalized electron 
Green's function and self-energy, while xio) n(g) are 
the fully renormalized spin-fluctuation Green's function 
and self-energy. The last term in (8) corresponds to the 
leading skeleton diagram of order 0{N) shown in Fig. 2. 
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FIG. 3; Illustrating functional derivative of the LW functional 
with respect to Green's functions. The cross indicates the line 
that is cut by functional differentiation. 



D. Eliashberg equations 

One of the important aspects of the LW functional ap- 
proach is that we can recover the self-consistent Eliash- 
berg equations for self-energies. Indeed, if we restrict 
ourselves to the leading 0{N) term in Y and use the 
stationarity of the free energy (8) with respect to self- 
energies (3), we get the following equations 
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+ 3g'TY,Gik + q)xiq)^ = 0, 
Il{q) + SNg^Tj^Gik + q)G{k) ) 



from which we deduce immediatly the expressions of the 
electronic self-energy and the collective spin polarization 

E(fc) = ?,g^TY,G{k + q)x{q), 
q 

Ii{q) = -3Ng^Tj2G{k + q)G{k). (9) 



These expressions can be also obtained by differ- 
entiating the leading order 0{N) contribution to the 
LW functional Y with respect to G and x respectively, 
according to equation (2). Diagramatically, this is 
equivalent to cutting one of the internal lines of the 
corresponding diagram, as shown in Fig. 3. 

Equations (9), with Dyson's equations, are nothing but 
the self-consistent Eliashberg equations for self-energies 
(see Fig. 1). 

Considering further terms in the LW functional Y 
amounts to studying deviations from the Eliashberg the- 
ory, in particular introducing vertex corrections as shown 
in Fig. 4, in the 1/iV expansion. 




FIG. 4; Spin-fluctuations self-energy generated from the con- 
tribution to the LW functional of order 0(1), shown in Fig. 2, 
by cutting one internal bosonic line. This can be obtained 
from The Eliashberg form of the spin-fluctuations polariza- 
tion shown in Fig. 1(c) by inserting a bosonic propagator to 
the vertex. 



E. Simplification of the Luttinger-Ward expression 

One can simplify further the full expression Eq.(8) of 
the free energy^. Indeed, from equations (9), one can 
notice that 

Y{G,x} = 3iV52T2^G(fc)x(g)G(fc-Hg) 

k,q 

= Arr^E(fc)G(fc) (10) 

fc 

= -T^n(g)x(9). (11) 

Thus, if we insert Eq. (11) into Eq.(8), the latter reduces 
to 



-ArT^[ln{-G-i(fc)} + E(A:)G(fc) 

k 

r^in{x-i(g)}, 



(12) 



while if we insert Eq. (10) we get the following expression 
for Eq.(8) 



Peff - - 



NTj2ln{-G-\k)} 
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+ 



TY^[]n{x-\q)}+Il{q)x{q)]. (13) 



Considering equation (12), we can show (See Appendix 
B 1) that the fermionic part reduces to a Fermi Uquid 

form FpL = — '^^P'' T"^ so that the final expression of the 
free energy for thermodynamics in the Eliashberg frame- 
work writes 

F,ff = FFL + TY^\n[x-\q)}- (14) 



F. Thermodynamics 

Performing the energy and momentum integrals in the 
EUashberg equations Eq. (9), one finds^^'^^ 



n(q,in) =7- 



(15) 



where 7 = Ng^xokp/iT^Vp). Considering only the linear- 
frequency Landau term in the spin susceptibility, this 
writes 

X-\ci,in) = Xo'{r' + |q- Ql^ + JS). (16) 

In particular, we have 2; = 2 for an antiferromagnetic 

QCP and z = 3 for a ferromagnetic one. 

Using this expression, the singular part of the free en- 
ergy (14) writes 

n\ 



/.(^^T) = T^ln{r^ + |,T+7^} 



1 

'27r3 



dv coth ^ 



2T 



dqif' tan~^ I 



IV I q'' 

(17) 



where q = q — Q is the shifted momentum near the wave 
vector Q. Performing the frequency and momentum in- 
tegrals in this equation, one finds the analytic expression 
of (17). Details of this evaluation for z = 2 arc given in 
appendix C. 

We would like to emphasize that this effective free en- 
ergy satisfies the following scaling relation 

fs{c\T) = 6-(''+^)^(^'6'/^^6^), (18) 

where fr{x,y) is an analytic regular function, u is the 
correlation-length exponent, z is the dynamical exponent 
and d is the space dimension. Inserting b = ^^"^ into the 
above scaling expression, we find 



ui,Ten- (19) 



Taking d = 3 and z = 2, the effective theory is beyond 
its upper critical dimension near the QCP, giving rise to 
V = 1/2, i.e., the mean-field behavior. Now, one can un- 
derstand thermodynamics near the HMM theory QCP 
based on this scaling free energy, derived from the effec- 
tive field theory in the Eliashberg framework. 



III. LUTTINGER-WARD FUNCTIONAL IN 
THE ELIASHBERG FRAMEWORK OF THE 
KONDO BREAKDOWN SCENARIO 

The HMM theoretical framework has been regarded as 
the standard model for quantum c;riticality in metals for 
a long time, although several heavy fermion compounds 
have been shown not to follow its predictions^^"^". An in- 
teresting alternative theory suggests that heavy fermion 
quantum transitions are selective Mott transitions of the 
/ nearly localized fermions^^"^^'^^ at which the Kondo 
effect breaks down. This scenario is supported by the 
presence of localized magnetic moments at the transition 
towards magnetism^'' and Fermi surface reconstruction 
at the QCPi*'!''. 

This problem has been tackled using the U(l) slave- 
boson representation of the Anderson lattice model^^'^'', 
with the introduction of a small dispersion for the /- 
electrons. A remarkable aspect of the theory is that 
the resulting QCP, at which an effective hybridization 
vanishes, is multi-scale. Indeed, because we have two 
kind of fermions in the model, i.e the conduction c- 
fermions and the /-spinous, there exist a Fermi surface 
mismatch q* = |/e^ — between Fermi momentum kp 
for spinous and kp for conduction electrons since fillings 
of spinous and electrons differ from each other. This mis- 
match gives rise to an energy gap E* for spinon-electron 
fluctuations that controls the dynamics of hybridization 
fluctuations. Although it depends on the value of q* , 
this energy scale is shown to vary from 0(10") mK to 
O(IO^) mK. When E < E* , holon fluctuations are un- 
damped, thus described by 2 = 2 dynamical exponent. 
On the other hand, when E > E* , holon fluctuations 
are dissipative since spinon-electron excitations are Lan- 
dau damped, thus described by z = 3 critical theory. 
Based on the z = 3 quantum criticality, recent studies 
have found quasi-linear electrical transport and logarith- 
mically divergent speciflc heat coefficient in d = 314,15,25^ 
and a divergent Griineisen ratio with an anomalous ex- 
ponent 0.7^^, consistent with experiments^^'^". 



A. U(l) slave-boson representation of the 
Anderson lattice model 

We start from the Anderson lattice model in the large- 
U limit 

^ = m 4a ((^T - - tij) Cj„ + ^ d\^{dr + ef)di^ 

i i 

+ V J2i4aCia + H.c.) + jJ2Si-Sj, (20) 

i {ij) 

where and di^^ are conduction electron with a chem- 
ical potential jj and localized electron with an energy 
level e/, respectively, tij the hopping term of the con- 
duction electron and V the hybridization between c- and 
d-electrons. The last spin-exchange term is generated by 
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a perturbative expansion to second order in t/U and is 
in competition with the hybridization term. 

In the U ^ oo hmit of (20), the strong correlations 
between the d-electrons show as a constraint of no double 
occupancy for the d-electron. This can be handled using 
the U(l) slave-boson representation 

d^a = bin,, (21) 

where bi and fi, are holon and spinon, associated with 
hybridization and spin fluctuations, respectively, obeying 
the local constraint 

^I^« + E/l/- = ^^' (22) 

a 

where S* = 1/2 is the value of spin and N is the number 
of fermion flavors with a — 1, N. 
One can then rewrite Eq. (20) into 

i 

+Vj2{hflc^,. + H.c.) + JY.{flx^,IJ. + H.c.) 

{ij) 

+NJ ^ Ix., + ^M^^ + - (23) 

The spin-exchange term for the localized orbital has been 
decomposed, using a fleld Xiji resulting in exchange hop- 
ping processes for the spinous. The local constraint (22) 
is taken into account by the introduction of a Lagrange 
multiplier A^. 

Performing the saddle-point approximation of bi — > b, 
Xij — > X) Etnd i\i — > A, one finds an orbital selective 
Mott transition as breakdown of Kondo effect at J « Tk, 
where a spin-liquid Mott insulator (b = 0) arises in J > 
Tk while a heavy-fermion Fermi liquid (6 =/= 0) results in 

Tk > J^^-^^. Here, Tk = Dexp(^j^^^ is the Kondo 

temperature, where pc ~ {2D)^^ is the density of states 
for conduction electrons with the half bandwidth D. 

Beyond the mean-field approximation, gauge fiuctua- 
tions corresponding to phase fiuctuations of the hopping 
parameter Xij = x^*"'^ should be introduced to express 
collective spin fiuctuations. It is more convenient to rep- 
resent the above effective Lagrangian as follows, perform- 
ing the continuum approximation, 

^ALM = Y I dv Cl{dr - tJic)Ca + ^^\drCa\^ 
+ fa{dr - - iar)fa + I (^i - iai)fa\^ 

2m/ 

+b*{dr - Mb - ^ar)b+ ^\{d, - ia,)b\' + 

+v{b*ej, + H.c.) 

+ ^f,..Uu + SN{iib + iar), (24) 



where g is an effective coupling constant between matter 
and gauge fields, and several quantities, such as fermion 
band masses and chemical potentials, are redefined as 
follows 

A -Mb, 
(2m,)-i-t, (2m/)-i = Jx, 
fi^ = fi + 2dt, - fif = €f + X-2Jdx- (25) 

In here, fermion bare bands and ej, for conduction elec- 
trons and spinous, respectively, are treated in the contin- 
uum approximation as follows 

el — — 2t(cos + cos ky + cos k^) 

« -2dt + t{kl + kl + kl), 

e{ = — 2 Jx(cos k^ + cos ky + cos k^) 

^~2Jdx + Jx{kl + kl + kl). (26) 

The band dispersion for hybridization can arise from high 
energy fiuctuations of conduction electrons and spinous. 
Actually, the band mass of holons is given by m^^ « 
NV'^Pcl2, where pc is the density of states for conduc- 
tion electrons^^'^^. Local self-interactions denoted by 
Ub can be introduced via non-universal short-distance- 
scale physics. One physical process for such interactions 
is four-point electron-spinon polarization (see Fig. 5), 

giving rise to Ub — uqj^ with uq ~ 0{\). Because 
such a local interaction term results from non-universal 
physics, one may consider that this term is introduced 
phenomenologically. 

/ i 9F i / 

FIG. 5: Four-point electron-spinon polarization for the 
holons. 

Maxwell dynamics for gauge fiuctuations appears from 
high energy fiuctuations of spinous and holons. 

We would like to develop the LW functional approach, 
including the Higgs or heavy-fermion phase. In this re- 
spect we write the holon field with its condensation part 
and fiuctuation contribution separately, 

b^B + b. (27) 

Then, the effective continuum Lagrangian is written as 
follows 

Calm j dr c*^{dr + :^9f ~ p,c)ca 
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-faiclr - ^J■f'- iar)fa- 



2m 



where 



+b* [dr - {fib - 2ubB^) - iar]b 



1 

2mb 



m~ia,)b\'' 



1 



+ y l&l' + Vib*c:U + H.c.) + VB{cU. + H.c.) 



( 

2mb 



2 



B^+(SN -B'^Y^ib + iar). 

(28) 



We see that the chemical potential for holon excitations 
is modified from fib to fib — 2ubB^. An important point 
is that gauge fluctuations become gapped when B ^ Q, 
which is due to the Anderson-Higgs mechanism. 



B. Luttinger-Ward functional in the Kondo 
breakdown scenario 



In the following, we demonstrate how thermodynam- 
ics can be extracted from the complicated effective field 
theory described by (28) , where two kinds of fermion ex- 
citations and two kinds of boson fluctuations are cou- 
pled with each other. The point is how to introduce all 
self-energy corrections self-consistently. As discussed be- 
fore, we construct the LW functional in the Eliashberg 
framework, allowing us to take all kinds of self-energy 
corrections self-consistently at least in the one-loop level. 

For simplicity, we start by ignoring gauge fluctuations 
corrections, considering only holon fluctuations. Gauge 
fluctuations are after that manipulated in the same way 
once their coupling with holons and spinous is known. 



1. Constructing a zero- order theory 

A subtle issue in deriving a LW expression for free 
energy is how to handle a non-vanishing condensation 
B ^ Q to describ the Higgs phase. A first step towards 
this derivation is to construct a " zero-order" theory tak- 
ing into account, in a proper way, the effect of the con- 
densation part B. 

Going to Fourier space, we can cast the action corre- 
sponding to the Lagrangian (28) into a mean-field part 
and holon fluctuations part : 



Smf ^-T^ [4fc5c ^{k)cak + flk9f ^{k, iuj)f„ 

k 

+VBTY, (/IfcC.fc + H.c.) + (^SN - B^)fib + Ub^ 

k 

Sfluc = ~T^b\d^'^{q)bq 

+VT^ E E + H.c) ' (29) 



9c ^(fc) =ii^ + fJ-c - 
9j^{k) ^iuj + fif - 



2mc ' 

k2 



2m 



/ 



d-\q) = in + fib - 2ubB^ - 



2mb ' 



(30) 



The interaction term gives raise to two kind of vertices, 
shown in the following flgures 
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FIG. 6: Vertices due to the interaction between fermions and 
liolons. Here, a line stands for the spinon propagator, the 
dashed line for the electron propagator and the wavy line for 
the holon propagator. 

Whereas it is justifled to follow the same strategy as for 
the HMM framework (See Sec. II C), i.e. use the cumulant 
expansion to the second order, for the fluctuations part 
of the interaction term, it is not the case for the conden- 
sation part. However, the latter can be considered as a 
renormalization of the propagators gc and gf and is thus 
included in a new zero-order theory whose bare action is 

Smf- 

Indeed, we can write 



where 



(^SN-B^)fib + Ub 
9 c' ~VB 



(31) 



-VB 



3f 



This gives the renormalized matrix Green's function for 
the fermions 



Go 



'^cc '^cf 



where 



9c 

gc'gf - {VBf 
9? 



gc'gj' - {VBY 



k,a q^O 



^0 _ ^0 _ 



VB 



gc'gj' ~ {VBY 



(32) 
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The condensation renormalizes thus the propagators 
for the f — f, c—c and / — c channels. In fact, this equals 
summing the infinite serie of the cumulant expansion due 
the condensation part of the interaction term (see Fig. 7). 



Y[G,x] 



^^^^ 



+ 




0{N) 



0{l) 



c f c 



c f c f c 



FIG. 7: Propagators of the zero-order theory, where the effect 
of the condensation B is totally taken into account. 



2. Derivation of the Luttinger- Ward functional 

Once we have properly handled the condensation part 
of the interaction, we can follow the strategy of Sec. II C. 

The interaction term due to hybridization fluctuations 
writes 



V -' » I, „ ,, „ /n 



(33) 



where we have extended the model to N identical species 
of fcrmions by adding a channel index i' £ [1, ^] to the 
fermionic operators and the 1 / \/N factor ensures a well- 
defined large-N limit. 

Considering only the leading 0{N) contribution to the 
LW functional shown in Fig. 8, and according to the gen- 
eral formula Eq.(l), we get the following expression for 
the free energy 



Ub- 



(34) 



with 



^TTr [In (-0(71 + 5])+ EG] 



Fb - T Tr [In {-d, 



Yb = ^2NV^ Y.k E,^o Db{q)GUk)Gff{k + q) 



(35) 



In (35), E 



are 



5- l^f ) and G = ( 

the self-energy and full Green's matrices, respectively, of 
the fermions. They are related by the Dyson's equation 



G- 



Gn 



The same equation holds for holons 



= d-' - Ub. 



(36) 



(37) 



FIG. 8: Leading skeleton diagrams participating to the LW 
functional Yb for the Kondo Breakdown model with depen- 
dence on 1/N. The fermionic propagators correspond to 
Gff (straight line) and Gcc (dashed line) Green's functions 
of the zero order theory including the condensation part of 
the interaction. Both fermionic and bosonic propagators are 
fully dressed. The first diagram contains one fermionic loop 
carrying spin and channel quantum numbers and one pair 
of vertices, each of order ©(l/x/iV), so that it is of order 
0{N^ /N) = 0{N). The second diagram contains four loops 
and four pairs of vertices, so it is of order 0{N'^ /N"^) = 0(1). 



One can manipulate gauge fiuctuations in the same 
way as the above, where the gauge-coupling action is 
given by 



k,q 
k,q',q 



2m 

k,q' .q 

whose vertices are shown in Fig-9. 



(b) 



(38) 





(c; 




(d) 




FIG. 9: Vertices due to the interaction of the gauge fields 
with holons and spinous. Gauge propagator is represented by 
a zigzag line. 
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3. Introduction of gauge fluctuations 

Following exactly the same procedure for hybridization 
fluctuations, one finds the following additional terms in 
Eq.(34) 

Fa = T Tr [In (-d-i + n,) + n^D,] 



F{q,k)Gff{k)Da{q)Gff{k + q) 



B{q, q')Db{q)Da{q')Dk{q + q'), (39) 



q,q' 



where 



"6 



(b) 



FIG. 10: First order skeleton diagrams corresponding to Ya- 
Fermionic and bosonic propagators are luUy dressed. 



F{k,q) and B{q,q') are the current-gauge bare vertices 
for spinous and holons given respectively by 



B{k,q)^^-j2vt{ 



q^ ' 



q^ J J"' ' ~ 7n, ' 



Ya corresponds to the contribution of the leading skele- 
ton diagrams, due to interactions with the gauge field, 
constructed with the fully dressed propagators of the 
spinous, the holons and the gauge fields (see Fig- 10). 



S,,(fc) = 2V''TY,Db{<l)Gn{k + q) 

= TYF{k,q)Gff{k + q)Da{q) 
q 

+2V^TYGccik - q)D,{q) 
q 

Ti,{q) = ni + ui' 

= TYB{q,q')Da{q')D,{q + q') 

q' 

+NV^TYGff{k + q)GUk) 

k 

iia{q) ^ n{ + n^ 

NT 



T 



YF{k,q)Gff{k)Gff{k + q) 

k 

YBiq',q)Dbiq)Db{q + q') (40) 



We see that the gauge field induces an additional part 
in the self-energies of the spinous and the holons, which 
we notice E^j and 11^ respectively (See Fig-11). 



5G 



ff 



Gff 
Gff 



Gff a 




G 



ff 



SYcL 




FIG. 11: Illustrating functional derivative of the gauge part 
Ya of LW functional with respect to Green's functions. The 
cross indicates the line that is cut by functional differentia- 
tion. Additional parts for the f-f channel and the holon self- 
energies are generated due to interactions with the gauge field. 



C. Eliashberg equations 

As for the HMM model, we can show that Eliash- 
berg equations can be derived from the LW functional 
approach. Indeed, restricting ourselves to the leading 
0{N) terms ofY = Yb + Ya shown in Fig-8 and Fig- 10, 
one can derive, in the same manner as in the HMM case, 
the following expressions for the self-energies 



Equations (40) are nothing but the Eliashberg equa- 
tions for the KB model studied in^^ . It is found that this 
approximation becomes exact in the limit N oo, pro- 
vided we take into account the Fermi surface curvature^*^. 



D. Simplification of the Luttinger-Ward expression 

We can notice that 
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k q^Q 

NTY,^cc{k)G,, + Y.)f{k)Gff (41) 

k 

-Tj^Kmiq), (42) 



and that 



Y„. = - 



^ F{q,k)Gff{k)D4q)Gff{k + q) 



y2 



2 ^S(9,g')A(5)^a(9')A(9 + 9') 
= 7VT5^E^/(fc)G/Kfc)-Enb(a)A(9) (43) 

k q 

= NTj2^f{k)Gff{k)-J2KiQ)Da{Q) (44) 

k q 

Hence, we can simplify the expression of free energy 
for the KB model by introducing either equations (41-44) 
into the expression of the free energy. As an exemple, if 
we use Eq. (42) and (44) we get the following expression 
for free energy 

Flw = FF + Fb + F,+ (^SN - (/.6 + iar) + u^^, 
with 

FiT = -T Tr [In (-G-^) + S^G] 
Ffe = T Tr [In {-D^') + A] 

F„ = T'IV[ln(-D-i)] (45) 

In this expression, the fermionic contribution reduces 
to a Fermi liquid form (see Appendix B 2) 



wNp. 



LW 



6 6 

where p± is the density of states of the upper(lower) hy- 
bridized band. The singular part of the free energy is 
given solely by the bosonic sector : Fg = Fa + F^. 

We can make a further simplification considering that 
n^'^ > n^, in which case the holon part is given by 

F, = TY,ln{-d^Hq) + nf%q)). 



E. Thermodynamics 

Performing the energy and momentum integrals in the 
Eliashberg equations [Eqs. (40)], we find Landau damp- 
ing expressions for the holon and the gauge polarizations 



U(%q,ifl) 



2m}) q ' 



n{(g,ifi)+n^(g,ifi) = ^M, (46) 



where 



76 



27r 7a iVTT 7r/d 



with nib = NWfT^ and fd = J (affrf-i) i with a UV 
cut-offi^ 

The singular contribution for the free energy is 



+TEln(g2 + 7,M + A,(M6,T)) 
+T^ln(g2 + 7„M+A„(/.b,T)), (47) 



where the condensation part fcif^'b^T), holon mass 
A(,(/i(,,T), and gauge-boson mass Aa{iJ,b,T) are given by 



Ab{p.b,T) = -2mb{p,b)[fj,b - 2ubB'^{iJ,b,T)], 

A^{^ib,T) = ^^B\^ib,T), 

rubilJ'b) 



(48) 



respectively. Note that the holon band mass depends 

on the effective chemical potential since it is given by 
the electron density of states. The coefficient in the 

gauge-boson mass is given by ^^I'/^^^o^ « 0{l)(j^'j , 
approximately. 

An important remark is that we can determine self- 
consistently the condensation value B by the condition 

= (See Appendix E). Beside the part obtained 
at the mean-field level, there are contributions due to 
hybridization and gauge fluctuation corrections. The 
Eliashberg framework allows then to reflne the value of 
the condensation B. 

An explicit analytic expression for the singular part of 
the free energy is obtained after integration on frequen- 
cies and momenta in Eq.(47) The details of this evalua- 
tion and its results are given in Appendix D. 



F. Scaling of the free energy near the Kondo 
breakdown quantum critical point 

Once we have the analytic expression for the singular 
part of the free energy, we can deduce its scaling expres- 
sions. As shown previously, this part of the free energy 
results from collective boson excitations associated with 
hybridization and gauge fluctuations. For each of these 
bosonic excitations, we associate a length scale for such 
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boson excitations, and the scaling form of the free energy 
near the Kondo Breakdown QCP reads 

+b-^''^^^^fa{\fbl/'^\Thl'). (49) 

fb{a){x,y) is an analytic regular function for hybridiza- 
tion (gauge) fluctuations and d is the space dimension. 
^6 = A^^''^ is the correlation length for holons, and 

\a = Aa is the one for gauge bosons. In particular, 
Ad may be considered as the penetration depth in the 
superconductor, bb and ba are scaling parameters for hy- 
bridization and gauge fluctuations, respectively, ly^a) is 
the correlation- length exponent of holons (gauge bosons), 
and Zh(a) is the dynamical exponent of holons (gauge 
bosons). Here, = and = Za are obtained. 

Although two kinds of length scales are introduced, 
both scales diverge at the same parameter point, V = Vc 
because they are related with each other via Anderson- 
Higgs mechanism. In addition, we note that this expres- 
sion is applicable near the Kondo breakdown QCP, ap- 
proaching from the heavy-fermion side because we have 
considered properly and in a self-consistent way the effect 
of a flnite condensation. 

Inserting bb = and ba = into the above scaling 
expression, we find 

/,(c^A-^T) 

= 4"''^''^'Vfc(i,Ter^) + A-2-('*+^)/a(i,TAri5o) 

Taking d = 3 and z = 3, the scaling free energy is written 
as 

/,(6, A„,T) = i^^fb{l,Teb) + A-Va(l,TAi),(51) 

where v = 1/2 coincides with the mean-fleld value owing 
to the upper critical dimensionality. Actually, the ana- 
lytic expression of the free energy turns out to follow this 
scaling relation. This is our main result, derived from the 
microscopic model based on the LW functional approach 
in the Eliashberg framework. Now, one can understand 
thermodynamics near the Kondo breakdown QCP based 
on this scaling free energy. 

IV. DISCUSSION AND SUMMARY 

In this study we derived the scaling of free energy from 
a microscopic model for two models of quantum criti- 
cality : the standard theoretical framework called the 
Hertz-Moriya-Millis (HMM) theory and the strong cou- 
pling approach corresponding to the gauge theory. Fluc- 
tuation corrections are taken into account systematically 
in the Luttinger-Ward functional approach. The Eliash- 
berg framework allows to use the proper level of approx- 
imation to get, self-consistently, the correct scaling for 
thermodynamics near the quantum critical point (QCP). 
We have shown that the singular part of the free energy 



for both models is due to the collective bosonic excita- 
tions, whereas the fermionic excitations give a Fermi Liq- 
uid contribution. 

For the HMM theory, there exists one length scale 
associated with the corresponding symmetry breaking, 
here spin-density-wave (SDW) instability. This fact al- 
lows us to construct the scaling free energy as a function 
of the spin-spin correlation length and temperature for 
the SDW quantum transition. We derived the correct 
scaling expression using the Luttinger-Ward functional 
approach in the Eliashberg framework. 

For the gauge theory, there are additional collective 
excitations. These have nothing to do with the phase 
transition directly although they are affected by it. Such 
collective modes turn out to be gauge fluctuations cor- 
responding to collective spin fluctuations in our context. 
An additional length scale, associated with gauge fluctu- 
ations, can appear. Indeed, considering that the Kondo 
breakdown transition is driven by condensation of holons, 
corresponding to the formation of an effective hybridiza- 
tion, the structure of the theory gives rise to massive 
gauge fluctuations via the Anderson-Higgs mechanism. 
This is the physical reason why the second length scale 
appears in the gauge theory. 

Because the two kinds of length scales, correlation 
length of hybridization fluctuations and penetration 
depth of gauge fluctuations, are deeply related via the 
Anderson-Higgs mechanism, they diverge at the Kondo 
breakdown QCP simultaneously. However, the presence 
of the additional length scale leads to a different scaling 
expression for the thermodynamic potential, compared 
with the HMM theory. We derived such a scaling expres- 
sion using the Luttinger-Ward functional approach in the 
Eliashberg framework. In the Eliashberg approximation, 
we showed that the scaling expression of the free energy 
has two contributions corresponding to each length scale, 
where each part contains only one length scale. 

In this paper we ignored vertex corrections, some- 
times justified but not always^°. Our path integral 
derivation of the Luttinger-Ward functional gives a 
chance to extend the Eliashberg framework, allowing 
vertex corrections. This can be achieved by going to 
higher orders of the cumulant expansion. In particular, 
if we do the same job up to the fourth order, we 
expect that vertex corrections will appear, satisfying 
the Bethe-Salpeter equation for vertices^"^. It is an 
important future direction to see how introduction of 
vertex corrections changes the scaling expression of the 
Eliashberg approximation. 

This work is supported by the French National Grant 
ANR26ECCEZZZ. 
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APPENDIX A: DERIVATION OF THE 
LUTTINGER-WARD FUNCTIONAL UP TO 
SECOND ORDER IN THE INTERACTION 

The LW functional can be derived thoroughly using a 
cumulant expansion to the second order in the interaction 
term. Indeed, this term induces the following corrections 
to the bare action 

k,k' q,q' 
k,k' q.q' 

(Al) 

The fermionic and bosonic propagators are introduced 



'9 -9 



' 2 



as 



on nn 



k + q,a 

FIG. 12: The one loop closed diagrams for free energy in the 
SF model. The fermionic and bosonic propagators are bare 



In Fig-12, the diagram (b) will generate a static and 
unifrom part in the self-energy and can thus be consid- 
ered as a renormalization of the electron chemical poten- 
tial. 



APPENDIX B: FERMIONIC CONTRIBUTION 
TO FREE ENERGY 

1. Spin-fermion model 

The momentum dependence of the electron self-energy 
is shown to be regular^^, we then replace the momen- 
tum with the Fermi momentum kp : = S(ki?,iw). 
Then, the electron contribution to the free energy of the 
SF model is given by 

F,i =-iVr^[ln{-Goi(A:,ia.) + SM} 



+i:{iuj)G{k,iuj) 



X{q)<>qq' 

while the corresponding self-energies are 

-g \Va'k'+q'Vi3fk'^al3'^-^q^a'f3''^~q' 
U{q)Sqq, = -g^y^ik + qT^liMlk' + q'^^'P'^P'k')^ 

The two last sums in (Al) corresponds to the two 
diagrams shown in Fig-12 where the fermionic and 
bosonic propagators are bare. 

The prescription to get the LW functional at this order 
is to dress these propagators^ in the corresponding skele- 
ton diagram, resulting in the following expression for the 
LW functional 

Y= N^4T^Y.k^^G{k)x{q)G{k + q) 

+ N^'4T^x{^)i:uG{k)Y.k'G{k'). (A2) 



= -NT 



du duln\^-Go\k,iuj) + uY.{iuj)^ 



7VT^ln{-Go-i(fc,zc.)} 



iuj^k 



iuj.k 



(Bl) 



Here, the last term corresponds to the free fermion part 
giving raise to the Fermi liquid form of the free energy 
for electrons 



-NT^ln\^-Go\k,iuj)^ 



iuj.k 



= -n'^T^, (B2) 



where is the density of states at the Fermi level. The 
first two terms of (Bl) are shown to be vanishingly small 
in the low energy limit. Indeed, we have 



6F,i = -NT 



du duln^-Ga\k,iLo) + uY.{iuj)^ 



S(iti)) 



Gq {k, iLj) — E(iaj) 



iu^k 



du 



GQ^{k, iLo) — uE(iw) 
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Gq {k, iuj) — S(?w) 



£-^Go-^(fc,za;)-SMio 



(l-») 
(B3) 



Now, we can switch from the integration over momentum 
to that over energy as follows 



Pf I de... 

-A 



where A is an energy cut-ofT. 
Integrating over e, we find that 



= -NTpF dw^^Szu; In, ). 

J /iw — MS(ici;) — A\ 1 
via; - uT,{iijj) + A/ J 



(B4) 



For |uS(ia;)| ^ A, the two last terms cancels and 5Ff,i 
vanishes. The electronic part of the free energy (12) of 
the SF model reduces then to the Fermi liquid contribu- 
tion 

= -n'^T^ (B5) 



2. Kondo breakdown theory 

The fermionic sector in the KB model factorizes into 
an upper (+) and a lower (-) band whose dispersions are 
given by 



ek + e£±A/(6k-e^)2 + 4F2B2 



The free energy for each band has a similar expression 
to Eq. (B5) and reduces to a Fermi Uquid form 

= -^T^ 

6 

where p± is the density of states of the upper (lower) 
band at the Fermi level given by 



P± = Po 



k± 



\E±=0 



APPENDIX C: MOMENTUM AND FREQUENCY 
INTEGRAL FOR THE FREE ENERGY OF THE 
SPIN-FERMION MODEL 



Introducing f{x) = tan ^ x, we have the following lim- 



its 



fix < 1) « X, fix > 1) 



(CI) 



Then, the free energy expression of Eq. (17) can be cast 
according to 



fsi^, T) ^-2^ C dy coth(^) ^ p dq 



+ g2 



, (C2) 



where qr 



^-2. Then 



/«(C,r) = -^^ rf^coth( 



2t) 



1 (2T)2 



3/2- 



27r3 



K+c'friTe), 



(C3) 



where Ag is a momentum cutoff and 



frire) = 



27r3 



212 



'2r$2 



- 1 



+ 



f- 

V2 



tan 



/2TC2 



We see that the singular part of the free energy follows 
the scaling relation shown in Eqs. (18) and (19). 



APPENDIX D: MOMENTUM AND FREQUENCY 
INTEGRAL FOR THE FREE ENERGY IN THE 
KONDO BREAKDOWN SCENARIO 

Let's consider the spectral representation of Eq. (47) 
fsilJir,T) = /c(Mr,T) - '^'^coth^— j dqq^ 

Considering the approximation for tan~^ x, the holon 
part is cast, as previously, into two parts in the momen- 
tum integral 



dv cothf I 
\2T) 



d^^coth(^) / dqq 



2 76«^ 



27r3 



dv coth 



f- 



2TJ.I, '^^^'i' 
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where Qr is a characteristic momentum determined by 
the equation 



76 



l^q^ + AbQr - -fbv = 0. 



The solution of the latter is given by 

(2/3)V3Ab 



+ 



/ , \ 1/3 

(976Z. + V12A3 + 81(76Z.)2j 



1/3 



21/332/3 



(D2) 



which is definitely positive. 
Then 



1 



dx 



2TJ Jg2 ""x + Ab 



127r2 



du coth 



A_^'^,,eoth(^)(-A,g. + 7..), (D3) 



where the momentum cutofi^ Ag is taken much larger than 
the holon mass, i.e., ^ Afe. 

The frequency integral can be performed approxi- 
mately, given by 

A(M.,T) = -^^°°d.coth(^)7..1n(^) 
-lii rf^coth(^)(-A,g. + 7.^') 

-1^1/ '^^V+L '^'^}(-^''«'- + ^^'^) 

1 /"'^ , 2T /A|g^AM^N 
/ di/ — 761^ In 1 



127r2 7o J/ 



'A2g,[A;„2T] 



) 



(D4) 



where 

9r[A6,2T] 



+ 



(2/3)V3A6 

(9[27bT] + Vl2A3 + 81(27fcT)2) 

(9[276T] + Vl2Ag + 81(276T)2) 



1/3 



1/3 



21/332/3 



(D5) 



For the gauge- fluctuation part, exactly the same pro- 
cedure is performed, and the result holds provided that 
the subscript h is replaced with a. 



APPENDIX E: SELF-CONSISTENT EQUATION 
FOR B 



Minimizing the free energy Eq.45 with respect to B, 
we get the following expression 



= Hmf{B) 



SrubUbB 



q2 + j,\f-2mb[l^b-2ubB^] 



+T 



y [— 



" y rrib 



)%^ + 7a^ + ^B2' 



(El) 



where Hmf{B) = determines the mean-field value 
of the condensation B. Fluctuations of the holon and 
the gauge fields result in additional terms in the self- 
consistent equation for B. 

The holon part is evaluated as follows 



1 







Ib^ + Ab 



d^q 
(27r)3 



ib^/q 



(g2 + Af,)2 + (76Z.)2/g2 



1 

2^ 
1 

4^ 
1 

47r3 
1 



00 



ifl — u 



dq 



du coth 



dv coth 



(2t) 

( 



dx 



q\q^ + Afc)2 -t- (7(,i.)2 

^hVX 

X{X + Afc)2 + {^^vY 
^bVX 



J dv coth ^ 



21"^ >/Maa;[A(,,(7(,!^)2/3] 



dx- 



47r3 Jo " \2TJ Max[Ab, (76i/)2/3] 

1 (276^)^ 
477376 Max[Ab, (276T)2/3] ' 



(E2) 



where the Max function is defined as Max[A,B] = A 
when A> B. 

The gauge part is evaluated in the same way and the 
analytic expression for the self-consistent equation of B 
writes 



= Hmf{B) 
rubUb 



+ 



B(27bT)2 



77375 Max[-2mb{iib - 2ubB^), (27bT)2/3] 
mg/mb g(27aT)2 
47r37„ Max[^e2, (27„r)2/3] 



(E3) 
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